New integrable model of correlated electrons with off-diagonal long-range order from 

so(5) symmetry 
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We present a new integrable model for correlated electrons which is based on a so(5) symme- 
try. By using an ^-pairing realization we construct eigenstates of the Hamiltonian with off-diagonal 
long-range order. It is also shown that these states lie in the ground state sector. We exactly solve 
the model on a one-dimensional lattice by the Bethe ansatz. 
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The study of models of correlated electrons is a significant tool in the theory of condensed matter physics. On a one 
dimensional lattice there are several known models which are exactly solvable by Bethe ansatz methods. The most 
famous of these is the Hubbard model whose solution was obtained by Lieb and Wu Q . Another well known example 
is the t — J model, the strong-coupling limit of the Hubbard model, which was in fact shown to be integrable at 
the supersymmetric point || through use of the Quantum Inverse Scattering Method (QISM) |3| . In this formalism 
the Hamiltonian of the model is derived from a solution of the Yang-Baxter equation, hereafter referred to as an 
i?-matrix, which provides a systematic method to obtain higher order conservation laws that guarantee integrability. 
An important aspect of the integrable coupling of the t — J model is that the i?-matrix is invariant with respect to 
\Q • the Lie superalgebra gl(2\l). For the case of the Hubbard model the symmetry algebra has been identified as so(4) 

os: §. 

A further important integrable correlated electron model was proposed and solved through the algebraic Bethe 
ansatz method by Essler et. al. ||. This model generalizes the Hubbard model with the addition of correlated 
hopping and pair hopping terms and is constructed from an i?-matrix invariant with respect to the Lie superalgebra 
gl(2\2). Another direction of generalization was given by Bracken et. al. || using the i?-matrix obtained from the one- 
parameter family of four dimensional representations of gl(2\l). The resulting model, known as the supersymmetric 
U model has also been solved and analyzed by Bethe ansatz techniques . In all the above examples the underlying 
symmetry has crucial consequences for the multiplet structure of the models providing insight into the ground state 
I i and elementary excitations. 

Recently it has been proposed that the antiferromagnetic and superconducting phases of high-T c cuprate compounds 
are unified by an approximate so(5) symmetry 0. Considerable support for this proposal came from numerical 
investigations in models for high-T c materials. In particular, it was shown that the low-energy excitations can be 
classified in terms of an so(5) symmetry multiplet structure Q. Subsequently extended Hubbard models related with 
an so(5) symmetry have been introduced and analysed in detail Jlo| . 

To our knowledge, no integrable correlated electron model associated with an so(5) symmetry has been proposed nor 
exactly solved. In this paper we construct such a correlated electron model which is exactly solved in one dimension by 
the Bethe ansatz. The integrability of our Hamiltonian comes from the fact that it is derived from an so(5) invariant 
i?-matrix which satisfies the quantum Yang-Baxter equation. Eigenstates of this Hamiltonian exhibiting off-diagonal 
long-range order (ODLRO) can be constructed through an ^-pairing mechanism. We also argue that these states lie 
in the ground state sector, which is a prerequisite for superconductivity. 
The Hamiltonian of this model is given by 
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H = J2h i , i+1 + h L , 1 + nN + BS z (1) 

where 



»=i 



m^c^ + h.c.) - 4(S*S* + SfS] + 3S?S]) 
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3(th + rij) - 3n l% - - 4[(l/4(n, - n,) 2 - (Sf - S*) 2 ] 2 (2) 
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Above Ci<j, cj CT are annihilation and creation operators for electrons of spin a, the Si's spin matrices and the Uf-f , n^'s 
occupation numbers of electrons at lattice site i. The number of lattice sites is L, S z — Yli=i &i * s the magnetization 
and N — X^i( n «T + n ii) 1S the number of electrons. This Hamiltonian exhibits correlated electron hoppings, pair 
hoppings, XXZ type interaction, chemical potential, nearest neighbouring Coulomb interaction, and the last term 
characterizes interactions favouring antiferromagnetism. The energy levels of the model are 

where the Uj's are solutions of the Bcthe ansatz equations 

(Ui + l/2\ L , % -i— i- Ui — U« I + 1 -T— r Ui — Uk — 1 
, = TT — 3 -^—Y\- ,i = l,...Mu 

i=n^-jn^i|,^i,...M 2 , ( 4 ) 

where Mi = 2L — N and M 2 = L — iV-p. Integrability of this model will be established through the QISM. The 
energy eigenvalues as well as the Bethe ansatz equations are obtained through the analytic Bethe ansatz jnj. The 
key ingredient to both of these methods is the following i?-matrix Jl2j 

R(u) = ( u ( u ~ 3)< ® e j + ( 3 " u ) e ) ® 4 + (-iy +3 u e) ® el j (5) 

which satisfies the Yang-Baxter equation 

R 12 (u - v)R 13 (u)R 23 {v) = R 23 (v)R 13 {u)R 12 {u - v). (6) 

Above the matrices e*- have entries (ej)f = SikSji, the indices range from 1 to 4 and i = 5 — i. This R matrix possesses 
the properties of 

- Unitarity 

R(u)R(—u) = (u 2 - l){u 2 - 9)1 <g> I (7) 

- Crossing-symmetry 

R u (u) = -A 1 R{3-u)A 1 (8) 

where t\ denotes transposition in the first space and A — ^-Y e \- 

The solution (^) is invariant with respect to the Lie algebra so(5) = sp(A) which has ten generators 

a i=e}-(-l) i +ld = -(-l) i +j<d (9) 

satisfying the commutation relations 

[a), af] = 6*a\ - 5}^ + (-l^fof - (-l) l+ ^a[. 

In order to build an electronic model we first need to put Z 2 grading in the i?-matrix. This is achieved by a 
redefinition of the matrix elements through 

RWii - (-i) W[jl+[fclM+ ™i?(^ (io) 

where we choose the parities to be 

[1] = [4]=0, [2] = [3] = 1 and [e*] = [i] + bl- 
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Such a matrix satisfies the Z2 graded Yang-Baxter where the multiplication of tensor products of matrices is governed 
by 

(o0i)(c8d) = (-l) lb]lc] ac<g)bd 
in eq. M). Following the QISM, we may construct the transfer matrix 



t(u) = str ^Rq L (u)Rol-i(u) ■ ■ ■ Ro2(u)Roi(u)j , (11) 

where stro is the supertrace over the zeroth space. From the Yang-Baxter algebra it follows that the transfer matrices 
t(u) form a commuting family and the associated Hamiltonian (FA) with /i — and B — can be obtained from 

H = r{uy l -t-t{u) 

du u=0 

where in view of the grading we have used the following identification 

|1) = |U>, |2) = |T), 13} |4), |4) = |0). 
In terms of the fermion operators, the so(5) generators (0) can be written as 



„t 



«i = -cj 



— 71 — 1 , (3) — 2S Z , Q>2 — ' — • ' — ( • ' ' > 

a\ = -c+, a{ = 2cW + , a\ = 2c+c-, al = 2S + , a\ = 2S~. (12) 
On the 2-fold tensor product space these generators act according to the co-product 

A(o}) = a) <g> / + (-1)™ ® a) for a) ^ a^a^a^af, 

A(aj) = a} (g) / + J® a) otherwise, (13) 

which extends to the i-fold tensor space co-associatively. Each of the local Hamiltonians ((U) are so(5) invariant. 

However due to the non-cocommutativity of the co-product the term breaks the so(5) symmetry of the global 
Hamiltonian (|l|). In spite of this, an so(4) symmetry is preserved comprising of an so(3) spin realization and an 
additional so(3) 77 pairing realization. For this reason we can add arbitrary chemical potential and magnetic field 
terms to the Hamiltonian which do not violate the integrability. 
The presence of the 77 pairing realization 



L L \ 

i=i j'=i j=i 

which can also be expressed in terms of the subalgebra generated by {a},a4,af}, allows a large number of states to 
be constructed exhibiting ODLRO |l3[| . Hereafter we treat the case [i = 0, B = 0. One can verify that if |0) = 
where |0) denotes the vacuum state. Thus the 2JV electron states 

|*Ar) = (vY |0> (15) 
are eigenstates of the global Hamiltonian with zero energy. These states are well known to possess ODLRO; that is 



lim - JV ' ^ J ' " J " =- 1- — (16) 

in the thermodynamic limit (Af,L — > 00, M/L fixed). Since the Hamiltonian is Hermitian the ground state energy 
satisfies 



E > LE, 







where Eq is the minimum energy of the two-site Hamiltonian. For this model we can determine that Eq = 0. It is 
thus concluded that the states dlq) lie in the ground state sector. 
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The energy levels (g) are determined from the eigenvalues of the transfer matrix (^Tj) which leads to a complicated 
expression that we will not give here. However, we mention that these eigenvalues are obtained through the analytic 
Bethe ansatz which exploits the properties of unitarity (]?]), crossing symmetry (^) and asymptotic behaviour of the 
i?-matrix. As usual, the Bethe ansatz equations are derived by the requirement that the eigenvalues are analytic 
functions. 

In conclusion, we have introduced a new integrable correlated electron model based on a so(5) symmetry. The 
model was exactly solved through the Bethe ansatz and shown to have ground states exhibiting ODLRO. 
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Desenvolvimento Cientifico e Tecnologico and Australian Research Council for financial support. J. L. and I. R. 
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